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' We derivate the Langevin and the Fokker-Planck equations for the radius of 0(3)-symmetric 

subcritical bubbles as a phenomenological model to treat thermal fluctuation. The effect of thermal 
noise on subcritical bubbles is examined. We find that the fluctuation-dissipation relation holds and 
that in the high temperature phase the system settles down rapidly to the thermal equilibrium state 
even if it was in a nonequilibrium state initially. We then estimate the typical size of subcritical 
bubbles as well as the amplitude of fluctuations on that scale. We also discuss their implication to 
ST ' the electroweak phase transition. 
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I. INTRODUCTION 



The dynamics of the electroweak phase transition is important in the light of electroweak baryogenesis 0. For 
the scenario to be successful the phase transition must be of first order with supercooling in order to attain a 
nonequilibrium state necessary for baryogenesis Though the effective potential of the Higgs field obtained by 
perturbation possesses a barrier between false and true vacua at the critical temperature || , it has not been certain if 
the phase transition accompanies supercooling. The purpose of this paper is to negatively confirm this. If the thermal 
fluctuation around the symmetric vacuum is too large, the phase transition proceeds without supercooling and no 
baryogenesis is expected. 

In the last five years, some attempts have been made by several different methods || - ||. For example, Gleiser and 
Kolb JjJ have discussed that the conventional picture of the first-order phase transition through nucleation of critical 
bubbles is applicable in the minimal standard model only in the case of relatively light Higgs mass, run < 70GeV, and 
that otherwise subcritical fluctuations plays a dominant role to realize emulsion of false and true vacua even above the 
critical temperature. In their analysis, however, it has been assumed that the typical scale of subcritical bubbles is 
given by the correlation length of the Higgs field. Meanwhile, as a way to understand the dynamics of quantum fields 
in a finite-temperature but nonequilibrium situation, Gleiser and Ramos [Q have derived the Langevin equation for a 
scalar field extending Morikawa's method jicj ■ The equation has also been numerically analyzed assigning white noise 
on a lattice, and it has been concluded that the sufficient phase mixing happen in any experimentally-allowed range 
of Higgs mass, where again the lattice spacing has been taken to be comparable to the correlation length [gj. Because 
the amplitudes of thermal fluctuation changes drastically depending on the length scale, it is important to determine 
physically the typical length scale of fluctuations which dominate the dynamics of the phase transition. As a first 
step along this line, in a previous paper |j| we have obtained the typical scale of 0(3)-symmetric subcritical bubbles 
by constructing the effective Hamiltonian for their radius and taking a thermal average. We thereby concluded that 
the electroweak phase transition is dominated by subcritical bubbles with any experimentally-allowed value of Higgs 
mass. Unfortunately, however, we cannot deny the fact that our approach was also too phenomenological to be viable 
from fundamental points of view. 

In the present paper, applying nonequilibrium statistical field theory we derive the Langevin and the Fokker-Planck 
equations for the Higgs field to clarify the detailed structure of the electroweak phase transition. Starting with the 



effective action we first derive these equations for generic field configurations. We then adopt a Gaussian ansatz 
for held conhguration to model a subcritical bubble and yield the effective Langevin and Fokker-Planck equations 
for its radius. This is a kind of variational approach which is adopted because, unlike the critical bubble, these 
subcritical bubbles do not constitute a solution of held equations but should merely be regarded as a model of held 
configurations. As a result our previous approach will be justihed under some conditions. Since we construct the 
Langevin equation in the perturbation theory, the effect from the environment is automatically taken into account. 
We need to calculate some loop corrections to obtain the friction term. The friction term will be estimated in the 
quasi-adiabatic approximation, while the noise term comes from the imaginary part of the effective action. The Fokker 
Planck equation can be derived directly from the Langevin equation and then we find that its static solution is of the 
form which we assumed as a probability amplitude in the previous paper. 

The rest of the present paper is organized as follows. In Sec. II, we review the non-equilibrium quantum held theory 
and discuss the statistical aspects of the theory. In Sec. Ill, we derive the Langevin and Fokker-Planck equations for 
the radius of subcritical bubbles. We examine the effect of dissipation and thermal noise on the electroweak phase 
transition. Section IV is devoted to summary and discussion. Throughout the paper we use the units c = U = 1. 



II. EFFECTIVE LAGRANGIAN OF HIGGS FIELDS IN THE THERMAL BATH 



A. Basics of the non-equilibrium quantum field theory 

In order to derive the Langevin equation and the Fokker-Planck equation for the subcritical bubbles, we use the 
non-equilibrium quantum held theory. This subsection is a brief review of its basic technique [ pT[ . 

The ordinary quantum held theory, which mainly deals with transition amplitudes in particle reactions, is not 
useful when we study statistical dynamics of macroscopic objects like bubbles. This is because we need the temporal 
evolution of some kind of classical order parameters with definite initial condition and not simply the transition 
amplitude of particle reactions with fixed initial and hnal conditions. The most appropriate extension of the held 
theory to deal with these issues is to generalize the time contour of integration to the closed form. More precisely, 
the time integration contour is generalized so that it runs from minus inhnity to plus infinity and then back to the 
minus inhnity again. This formalism, often called as in-in formalism of quantum held theory, yields various quantum 
averages of operators evaluated in the in-state without specifying out state. On the other hand the ordinary quantum 
held theory, often called as in-out formalism of quantum held theory, yields quantum averages of operators evaluated 
with an in-state at one end and an out-state at the other. 

The partition function in the in-in formalism for a real scalar held is defined to be 



Z[J] = Tr 

= Tr 
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where the suffix c in the integral means that the time integration contour runs from minus inhnity to plus inhnity and 
then back to the minus inhnity again. All the held quantity is defined on this closed time contour. In the above, X + 
represents a held component X which is restricted on the forward branch (—00 to +00) of the time contour and A_ 
stands for that restricted on the backward branch (+00 to —00). In the rest of this paper we often use the following 
notation: A a = X + — X- and Xq = (X + + A_)/2. The symbol T designates the operator ordering with respect to 
this closed time contour, and accordingly T + designates the ordinary time ordering and T_ the anti-time ordering. 
Here J means the external held. Although J + and J_are the same actually, we regard that they are different from 
each other for technical reasons and we only set J+ = J_ at the end of calculations. The symbol p represents the 
initial density matrix, and the held <f{x) is in the Heisenberg picture. 
In the interaction picture, this partition function becomes 
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where we have used the Wick theorem which holds not only in the vacuum state but also in the thermal state with 
p = exp(— H/T). In the latter the above propagator is, in the momentum representation, 



p 2 — m 2 -\-ie 



2irin(p)5(p 2 - m 2 ) -2m(0{p o ) + n(p))5(p 2 - m 2 ) 



-2m(6{-p )+n(p))5(p 2 -m 2 ) 
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(2.3) 



where n(p) is the thermal distribution function: n(p) — [exp(w(p)/T) — l]^ 1 and w = \/p 2 + m 2 . 

By the Legendre transformation of the partition function, we obtain the generalized effective action, or generating 
functional of the vertex functions: 
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immediately follows as in the in-out formalism and this form is often used as a generalized classical equation of motion 
for the variable <fr(x). However, this effective action has an imaginary part. For example if we parameterize the kernel 
of two-point part of T as [O 



rW(x,y) = 
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the imaginary part of I^ 2 -* is 
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We can rewrite this expression by introducing a real auxiliary field £(.t) JK 



exp( l r[0]) = J [^]P[e]exp 
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where, 



P[£\ = exp 



(2.9) 



is a normalizable positive definite statistical weight for the fields £(x). If we apply the variational principle for the 
exponent of the integrand of eq.(2.8), we obtain an equation of motion for (p c (x) as 
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This is a Langevin type stochastic differential equation with a nonlocal kernel A(x — x'). If the time scale of change 
in (f> c (x) is small compared with that of radiative corrections and <j) c (x) is nearly homogeneous in space, this term 
reduces to the familiar friction term: 



U4> c {x) + Vi s {4> c {x)) + r]<j> c (x) = 



(2.11) 



with 
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... dA(k) 

where V c s is the effective potentiai and A(k) is the Fourier transform of A. 
If we define the statistical average as 



= J [de]p[e]..., 



we then obtain 



{Z(x)ttf))e=B(x-J). 
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Thus we can construct consistent statistical field theory in the in-in formalism of quantum field theory. 

In the actual application of this formalism to the subcritical bubbles, we need to calculate the dissipative and 
diffusive kernels A and B at finite temperature. This is the subject of the next two subsections. 



B. The fluctuation-dissipation relation and the stationary distribution for Higgs fields 

For simplicity, we consider the following Lagrangian of a singlet Higgs field if. 



(2.15) 



Strictly, we must include interactions with Z, VK-bosons and all quarks for the detail of the electroweak phase 
transition. In particular, contributions of Z and W boson are crucial to yield the cubic term in the effective potential 
with one-loop corrections. However, the essence of non-equilibrium phase transition can be fully obtained by the 
above simple model as we will see soon. According to the calculation by Morikawa j[o| and Gleiser et al the 
effective action becomes, up to one-loop corrections, 

A 
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Here S^(x — x') is the thermal Green's function of the fermion defined by 
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where /(p) = (exp(|p|/T) + l)"^ 

We can rewrite this expression by introducing two real auxiliary fields £1 (x) and £2 (%) ■ 



exp(»T[ ¥ >])= / [dei][«a]fi[6]ft[6]exp iReT + i /(&cM+6<Ma) 



where 



P&] =exp 
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Applying the variational principle on eq. (2.23) as in the previous subsection, we obtain the Langevin equation: 

U<j) c {x) + V^(4> c (x)) + 2 [ d 3 x' I dt'^x^A^x - x') + <f> c (x) f d 3 x' f dt'(f> 2 c (x')A 2 (x - x') 

= 6(z)+0c(*)k(z). (2.25) 
Assuming that 4>{x) is nearly homogeneous and changes slowly in time, we obtain 

U(j> c {x) + Vl s {<t> c {x)) + r)x4> c {x) + m 4>l{x)$ c {x) = Zi(x) + M x )&( x ) (2-26) 

with 



and 



V2 = ^J dtj d 3 xA 2 {x, t)t, 



(Zi{x)Zi{a/)) £ = Bi(x-a/). 



(2.27) 
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For the contribution from only the self-coupling, Gleiser and Ramos j7| have obtained the result 772 = ^ m (^:) f° r 
the friction coefficient and 



(2.30) 



in the high temperature limit. However, for A ~ 10~ 2 (Higgs mass mu ~ 60GeV) the correlation time of the noise 
is A£ no ; sc ~ (Decay width) -1 ~ 15367r(A 2 T) _1 ~ 10 7 T _1 which is much larger than the typical scale ~ T _1 . Hence, 
the above approximation is not guaranteed in 

In our present analysis, on the other hand, the most dominant contribution comes from the Yukawa interaction 
with 1> / 3> A, so that terms proportional to A 2 are negligible and it suffices to consider the following Langevin 
equation: 



□^c + Vrfr(&)+»7i^c = 6i- 



(2.31) 



To obtain the friction term we must prepare the full propagator of fermion. Up to one-loop order, the renormalized 
Green's function becomes 



S F (p,t) 



2ld n 



(-Wft) + 7 ' p) /(-w p + iT f )e-^M 



7 °o; p e(t)+7-p)/(w p + ir / ) e ^l*l , 



(2.32) 



"Here we use a massless propagator as an approximation. In fact once the phase transition starts, the top quark has a space- 
time dependent mass due to the variation of the Higgs field. The maximum value of the top mass, however, remains as small as 
f<f>o ~ 50GeV at the critical temperature T c ~ 93GeV. We may therefore conclude that the result would not change significantly 
even if we used the more complicated propagator with nontrivial space-time dependence or simply a massive propagator. 
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where T / is the decay width given by 



8tt 

in the high temperature limit. The friction coefficient then becomes 



T f ( Wp ) = - 7 °Im(E) \ P 2 =0 ^^-T (2.33) 
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Here we are using the approximation that the external momentum is zero. In the case of Yukawa interaction with top 
quark, / <~ 1 and then /3Tf <~ 0.04. Thus, we can expand the last expression by j3T f and to the lowest order we find 

771 * 8^ff coBhOfc*) + 1 = 3 T (2 ' 35) 
Taking the same approximation used in the derivation of eq.(2.34), the autocorrelation of noise becomes 

■^S aP (k, t - i')^ Q (k, t' - t)\ <5 3 (x - x') 

= f i ^ 4[cosh ( ;^tIs( /3 r /) ]4 2co ^ r ^[ cosh (^) + 

- 2sin(/3w fc )sin(/3r / )sm(2w fc |t - t'|)}<5 3 (x - x'). (2.36) 
In the limit Tf » 1 (strong coupling limit / 3> 1), 

f 2 /■ d 3 k 1 r 
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- 2sin(/3w fc )sin(/3r / )sin(2w fc |t - t'\)^5(t - t')5 3 (x - x') 
I 2 f°° 



4ir 2 T f J cosh(/3w fe ) + 
~ 2r h T5{t-t')5 3, {yL-yi). (2.37) 

Properly speaking, one cannot take the above limit in the minimal standard model. But, we expect that the ap- 
proximation does not give drastic changes on the final result as long as we concentrate on a time scale larger than 

L f ■ 

Since the ordinary fluctuation-dissipation relation holds as above, we expect that the static solution is the canonical 
distribution. Introducing a new variable 

( 2 - 38 ) 



we rewrite eq. (2.31) as 

dp<j> c 



dt 

The Fokker-Planck equation for the above Langevin equation becomes 
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where 
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and J is the probability current density with the expression 

1 



J 



1 -^J\& W + T S: 



Here n is the Hamiltonian given by 
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From the expression of the probability current density we easily find a static solution of the system as 
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On the other hand, normalizable dynamical solutions may be expressed in the form | fl3| 
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The lowest eigen value Ao is of course zero, while the next lowest eigen value gives the time scale of relaxation to the 
thermal equilibrium. For its estimation we only have to consider the Hermitian part of (2.44) p3| , 



AW*W = / d 3 x2 m T 



2<5p 0c 2 (x) 



2T 



4T 
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The above equation is an infinite collection of harmonic oscillators with their ground-state energy subtracted. We 
thus find formally A^ — rji, that is, the time scale of relaxation is in general given by the inverse of the friction 
coefficient in the original Langevin equation. The relaxation time reads Ai re i ax = Re(Ai) -1 > A^ 



-l 



m c J ~ 10 2 GeV ~ 10 27 s for the electroweak Higgs fields. On the other hand the expansion time scale of the 
universe in this epoch is 



At 



expand 
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H 



10- I2 s. 



As At re i ax <C Ai oxpan d, the system is almost always in thermal equilibrium. 

By integrating over p^ c in (2.42) we obtain the probability distribution function P s t[4> c ] for 
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(2.48) 



where T4ff(0 c ) is the effective potential with loop corrections. Thus, we find the ordinary expression for the probability 
distributional function. 

We now obtained the basic tools for understanding the dynamics of the electroweak phase transition. In the next 
section, we estimate the amplitude the thermal fluctuation by adopting the Gaussian ansatz for the subcritical bubbles. 
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III. THE DYNAMICS OF SUBCRITICAL BUBBLES 



A. Derivation of the Langevin Equation for the Radius 



Having derived the equilibrium probability distribution function (2.42), we can calculate the characteristic spatial 
scale of thermal fluctuation by adopting, say, a Gaussian ansatz for the Higgs field as was done in the previous paper 
||. However, since we are interested in the kinematics of subcritical bubbles, we adopt a variational principle approach 
to analyze the effective action directly. We choose a trial functional of the form 



Do exp 



Rl(t) 



(3.1) 



to derive the effective Langevin equation for the radius R(t). The above ansatz is reasonable because its 0(3) 
symmetry helps to minimize the free energy and we are treating bubble with a thick wall at or above the critical 
temperature. Here (f>o should be identified with a local minimum of the effective potential which would become the 
global minimum below the critical temperature and is the most expected value apart from <f> = 0. We shall investigate 
the kinematics of subcritical bubbles with amplitude (f>o . The detail discussion for the above ansatz has been done in 
the previous paper ||. 

Thus the degree of freedom has reduced to one and we then insert eq. (3.1) into the effective action (2.16). 
First we calculate the real part. The result is 



Rc[r(i?+,i?_)] = 2tt 3/2 ^ / dt 
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where 



A(Kt-t') :=/ 2 



d 3 i 



(27T)'- 



-Re 



iSL,(p, t - t')S$ a (p -k,t'- t) 9(t - t') 



(3.2) 



(3.3) 



The imaginary part which generates the noise term becomes 

Im[T(R+, R.)] = ^J dtdt' J tSt (t)e-* fc3ii +W - Rt(t)e~^ R ^) (R 3 + (t')e-^ 2R ^ (*'> - Rt{t')e'^ k 
x B(k,t-0, 
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isZ p { Vl t-t>)s p F a {p-k,t> -t) 
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(3.4) 
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This can be rewritten with an auxiliary fields £(k, t): 

d 3 k 



exp{i x ilm\T(R + ,R-)]} = / [d£]P[£]exp 
where 

P[£]:= iVexp 
Thus the effective action for R± (t) becomes 



-- I ilh.it' 



(2tt) 3 / 2 
d 3 kd 3 k 



dt£(k,t)(R + (t) 3 e-i k2R lW -i?_(t) 3 e-5 fe2 «-(*))l, (3.6) 



(2tt) 3 
d 3 k 



£(k, t)B-\k, t - t')5 3 {k - k')£(k', t') 
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(3.7) 
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From SS c s/SRa\r & =o = 0, we obtain the effective equation for R c (t): 

cP R c 1 , dR c . 9 2 1 2 9 A (9 . „ 
+ 77TT(-^) + ^TT + + TT^O + ' ' -)Rc 



dt* 2R c y dt ' ^5i? c ^V" 6(VT° 

16 ^_| dt 'J _0_^( k ,t - f )(3i2e(t) - fc 2 ^(t))e-^ 2 ^Wi?3 ( ^ )e -^ 2 ^(t') 



+ 



15tt 3 / 2 ^, 

157t3/ 2 ^7 (2^)3/2 V 



3flc(t) - fc 2 i^(t))e-* fc ^ (t) ^(k,t). (3.9) 



Assuming that the time-dependence of R(t) is weak, the last term in left-hand-side of this equation is expanded as 



16y/2 f* f d 3 k 



+ i^^/./^"^ / ^(k,*-*')(3-fcX 2 (*))V^^W. (3.10) 

The first and second terms in the last expression give the loop correction on the effective potential and friction term 
Ff, respectively. Consequently Fj ~ r]R c , where r\ ~ r)i/4\/2. 
Next we turn to the noise term, which we denote as 

Ht) := I J^y ^(aflcW - ^(^"^^(M)- (3.11) 

By using the definition £(k, t) the self-correlation under the limit » 1 is given by 

(£(k, i)£(k', *')) = B(k, t - 0<S 3 (k - k') (3.12) 
and then the correlation of £(t) becomes 

(l(t)«*'» = 225^gRg l^/ ^ {X}( ' {x ' y > H W " ^(<>>-*"*'«"(3^(i') - ^(f))^-*^'''' 

In the above result one might suspect the i? c -dependence of the diffusion coefficient, but one shall understand by the 
later discussion that it is a quite reasonable result. 
Consequently the Langevin equation becomes 



where r\ := j^y^T and 



d^ R r 1 .dR r . 2 1 /m , „ dR r ~, . .„ ^ 

V + 2s:^) a + 5 + a(W = -our + ^ (3 - 14) 



a(T) = | m » + ^ + .., (3.15) 



B. The Fokker-Planck Equation for the Radius and its Static Solution 

In order to find the equilibrium state we transform Langevin equation cq. (3.14) to the Fokker-Planck equation. 
The Langevin equation can be rewritten as 

§ = Si - \i + - » p = + M *>- < 3 ' 16 > 
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where P c is the canonical momentum defined by 



Pr := 



<9L e ff _ 15tt 3 / 2 
8V2 



dRc. 
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(3.17) 



with M(T, R c ) :-- 



15tt 



3/2^2 R 

— - and L c ff is the Lagrangian of the deterministic part of the Langevin equation: 



8/2 



L eS (R c , V c ) = \mV 2 - 2 -M - X -Ma{T)R\. 



(3.18) 



Now we introduce the probability distribution function W(R C , P c ; t) of R c and P c . Then we obtain the following 
Fokker-Planck equation from the above Langevin equation: 



Pc 

M 



W 



d 
dP c 



p2 



2 M 



2R C M 5 R c 



Ma(T)R c + V P C 



W 



K,2~ d 2 
2 V 8P? 



W. 



The above equation can be expressed by using the effective Hamiltonian 



H c g(R c , P c ) :— P c Rc ~ L e g 



2M' 



as 



The probability current density becomes 



-M + -Ma(T)R 2 c 



dH c ft dH c ff 
+ Mr] 



dRc 



dP c 



W 



1 , «-9 . u 



w. 



(3.19) 



(3.20) 



(3.21) 



J = 



1 

-Mi] 



dHtf! TT/ , rp dW 

0P C 9P C 



and eq. (3.21) reduces to <9 t W / + V- J = 0, where V = {dn c ,dvj- Here we have used the relation between the diffusion 
and friction coefficients: 



^M 2 ij = TMr]. 



Then the static solution becomes 



W st {R,P) cxexp 



(3.22) 



(3.23) 



As in the case of Higgs fields the relaxation time scale is equal to ry -1 and is again much shorter than that of the 
expansion of the universe. Thus the distribution is always in the stationary state. 

Finally, we note the fluctuation-dissipation relation exactly holds. Let us recall the case of the Brownian motion of 
a massive particle. In such a system the relation between the diffusion and friction coefficients is 



(Diffusion coefficient) = (Friction coefficient) x (Temperature) x (Mass). 



(3.24) 



As the mass term corresponds to M in the present case, the above relation also holds. The i? c -dependence of W s t is 
a direct consequence of the three-dimensional volume effect. Further this dependence is physically reasonable. When 
one compares large bubbles and small bubbles, the destruction due to the thermal noise is not effective for large one. 



C. The Size of Subcritical Bubbles and the Amplitude of the Thermal Fluctuation 

Now we can estimate the amplitude of the thermal fluctuation. We first calculate the averaged radius of the 
subcritical bubble. Note that previously all authors have assumed the correlation length for the typical spatial scale 
of the thermal fluctuation. As we will see soon, the scale is a calculable quantity in high temperature phase. In this 
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subsection, we extend the above analyses to the electroweak phase transition. For this purpose we start from the 
following effective Hamiltonian corresponding to eq.(3.20): 



H eS (R c , P c ) = -^Pl + + ^Ma cw (T)Rl 



where 



and 



a C w(T) := -£>(T 2 - If ; 



3ET r 



8V2 
15V3 



10V2 



l-|^(T 2 -Tf)]. 



2A T '" T 9£ 2 T 2 
In the above expression, some coefficients can be fixed by electroweak particles: 



6( !^ )2 + 3(^) 2 + 6(^) 2 
a a a 



0.169, 



E = 



12tt 



6(^)3 + 3(^)3 

a a 



0.00965, 



\t = A — 



167r 2 L^s y V cr ' c B T 2; ^F yty a ' cfT 2 ' 



0.0350, 



(3.25) 
(3.26) 

(3.27) 

(3.28) 
(3.29) 
(3.30) 



and 



D 



64tt 2 



T? = ^/(m 2 H -8Ba 2 )/4D, 



6(^ + 3 (^)^l2(^ 
a a a 



-0.00456, 



(3.31) 



(3.32) 



where we used mw — 80.6GeV, mz = 91.2GeV, m f = 174GeV and a = 246GeV |L4j. Further we assumed m# = 
60GeV which is the minimum value experimentally allowed. 

Let us sketch the potential change with temperature. As the temperature decreases, the non-symmetric vacuum 
appears at T = Ti := T^/yl — 9E 2 /8XtD ~ 93.52GeV and then the hight of the two vacua coincides with each 
other at the critical temperature T = T c := T^/ ' \f\ — E 2 /XtD ~ 93.43GeV. Finally, the barrier between two vacua 
disappears at T = T%. 

As there exists the interaction with other particles in the present case, the expression of the friction and diffusion 
term obtained in the previous section should be modified accordingly. However, as long as the fluctuation dissipation 
relation holds, which has been confirmed in the presence of bosonic interaction in j?]] and fermionic interaction in the 
present paper, the distribution of Higgs fields is given by 



W st oc exp 



ffoff 

T 



(3.33) 



Because the system is always in the stationary state, the typical scale can be calculated by the ordinary canonical 
averaging. In the high temperature phase the thermal average of R is given by 



(Rc) 



JdP c dR c R c exp 
JdP c dR c exp[~^] 

J dR c R?J 2 exp [— — M £f™ R 2 } 
^dR c RV 2 exv[-W ~ 



$™dR c Rl /2 e^[-?§] 

J™dR c Rl /2 eM-W} 
2V2 T 



(3.34) 
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In the low temperature phase, a ew (T) is negative and accordingly the averaged radius diverges. However, this is not 
relevant. We are interested in whether the system has supercooling or not in the high temperature phase only. 

Having seen that R is the important scale of subcritical bubbles near the critical temperature, we shall now estimate 
the amplitude of the fluctuation of <j) around 4> = on this particular scale, adopting the trial function 



<j) = ^Aexp 



Then the free energy becomes 
r37r 3 / 2 



F( 



(R) 



r 3/2 



-_D(T 2 - t£)(rY 



4\/2 2y/2 
Thus the RMS amplitude of <f> at the symmetric vacuum is 



3%/3 



ET4>\{Rf 



r 3/2 



32 



A2 _ 



A<t> 2 A e- 



\ Jd<t> 



A<i~ 



16D(T 2 ~T 2 )T 2 



(3.35) 



(3.36) 



(3.37) 



Its temperature dependence is depicted in Fig.l. In the same way as in the case of (R), the above approximation 
result gives the lower bound for \J~^ 2 - At T — T c and T = Ti, one obtains numerically \f(jP(T — T c ) — 36.2GeV 



and 



Ti) — 28.2GeV, respectively. This exceeds the first reflection point 



ET 

At 



E 2 T 2 2D(T 2 -T 2 ) 



which implies that the perturbation theory breaks down and that one can no longer conclude that electroweak phase 
transition is of first order. 

The above argument implies that the phase transition is not accompanied by supercooling. The transition rate 
between two vacua is much larger than the expansion rate of the universe. In fact, it is given by 



77i#exp 



T 



(3.38) 



where F is the free energy of subcritical bubble configuration. As F ~ T during the high temperature phase M, 
r ~ 10 2 GeV ( ^> H ~ 10~ 13 GeV). Hence the fraction of two vacua are determined by the hight of the potential and 
therefore the fraction is almost the same at T = T c . This means that any supercooling does not happen and critical 
bubbles cannot borne. 



IV. SUMMARY AND DISCUSSION 



In the present paper we have examined the relaxation of the system from non-equilibrium states and found that 
the system settles down rapidly to the stationary distribution, which is simply the ordinary thermal distribution. For 
the Higgs fields it is 

P^cxexp^-ffl), (4.1) 

where T[4>\ is the free energy. For the radius of subcritical bubbles, 

W«{R,P)KeKp(-Z2&&y ( 4. 2) 

where H c g(R,P) is the effective Hamiltonian. First we estimated by using W s t(R, P) the mean radius of subcritical 
bubbles which determines the spatial scale of the thermal fluctuation. Next we estimated the field fluctuation by using 
P s t ■ We conclude that the electroweak phase transition is quite weakly of first order and therefore the standard 
baryogenesis does not work in minimal standard model, apart from the too small magnitude of CP violation. 

Writing down directly the Langevin equation for the radius we have a glimpse of the kinematics of subcritical 
bubbles. In the fluctuation-dissipation relation the dependence on radius is a natural result from the three volume 
effect. For larger subcritical bubble the effect of diffusion is smaller. 
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Finally we discuss validity of some approximations employed here. First in the derivation of the Langevin equation 
(2.31) we assumed the approximate homogeneity on the Higgs fields. In the derivation of the Langevin equation 
(3.14) for the radius we have used a similar approximation, which is correct if T > (R c ) . In fact (R c ) /T ~ 
(E/Xt) 2 ~ 0.1. Thus the approximation is justified. Second, we have taken the strong coupling limit to obtain the 
white noise. In practice it is colored and in general the diffusion term is written by a integral expression, which would 
be too complicated to treat analytically. We expect, however, as long as we focus on a time scale larger than TJ 1 the 
more elaborate analysis of the colored noise would yield the same result and that in this sense our simplified analysis 
suffices. We hope to consider this subject further in a future publication. Third, we did not work out the calculation 
of the loop correction by gauge bosons on the friction and diffusion coefficients. Although it is essential in order to 
obtain the cubic term in the effective potential, their contribution to the friction and diffusion coefficients is smaller 
than top quark. 
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Figure Captions 

1. Fig.l: The RMS of thermal fluctuations. The unit of the the vertical axis is GeV. 
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